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ABSTRACT 


A  differendai  equation  describing  the  energy  transfer  between  a  fluid  and  a  body 
moving  in  that  fluid  was  derived.  The  derivation,  based  upon  the  Coriolis  form  of  the 
NavieT'Stdces  Equadon,  contains  a  rigorous  viscous  correcdtm.  For  inviscid  i^al  cases, 
the  equadtm  demonstrates  that  the  rate  of  total  enthalpy  transfer  from  (or  to)  the  system  is 
a  function  of  the  transverse  component  of  the  pressure  gradient.  Therefcne,  for  practical 
turbomachinery  rottns,  the  derivative,  dp/dO  ,  can  never  vanish. 

On  integradon  of  the  differential  equations,  a  form  of  the  Euler  IXnbomachinery 
Equadtm  with  viscous  correction  is  derived.  The  resultant  form  contains  two  distinct 
work  rate  terms  for  the  axial  and  radial  componeu^s  of  the  flow.  The  fact  that  integratitm 
yields  a  result  which  approximates  the  classic  Euler  Turbomachinery  Equadon  consdtutes 
confirmadon  of  the  d  ^vadon. 

An  applkadon  of  the  equadon  to  an  ideal  infinite  linear  cylinder  with  bound  vcndc- 
ity  was  developed,  yielding  the  expected  known  result 
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NOMENCLATURE 
any  continuous  function 
radial  fraction  of  the  mass  flow  rate 
axial  fractitm  of  the  mass  flow  rate 
enthalpy  per  unit  mass  of  fluid 
total  endia4>y  per  unit  mass  of  fluid 
total  endialpy  of  a  flmd 

total  enthalpy  of  a  fluid  per  unit  length  of  blade 

relative  total  specific  rothalpy 

lift  per  unit  length  of  blade 

mass  flow  rate 

radial  mass  flow  rate 

axial  mass  flow  rate 

local  static  pressure 

total  pressure 

a  generalized  curvilinear  coordinate 
a  specific  heat  rate 

radial  coordinate  of  a  cylindrical  coordinate  system 

entropy  per  unit  mass  of  fluid 

time 

absolute  temperature 

velocity  of  the  blade  and  a  function  of  die  radius 
velocity 

the  integrated  viscous  term  of  Equation  (22) 
relative  fluid  velocity  in  a  moving  rotor  frame 
a  time-dependent  component  of  the  velocity  in  the  moving  frame 
axial  coordinate  of  a  cylindrical  coordinate  system 
the  circulation 

tangential  angle  coordinate  in  cylindrical  coordirates 
die  kinematic  viscosity 

the  stress  tensor  separated  from  the  thermodynamic  pressure 
local  fluid  density 

volume 

a  stream  function 
angular  velocity  vector 
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INTRODUCTION 


A  monument  of  turixtmachinery  technology,  the  Turbomachinery  Equation  of  Euler 
is  based  tq>on  thermodynamic  defmitions  of  work  and  Newton’s  Laws.  Since  the  Navier- 
Stokes  Eldons  and  Crocco’s  Equation  [1^2]  in  a  rotating  (moving)  frame  are  also  based 
upon  thermodynamics  and  Newton’s  Laws,  they  must  in  principle  contain  the  TUrboma- 
chinery  Equation  in  differential  form  and,  on  integration,  in  integral  form.  Application  of 
Coriolis’  transformation  to  the  viscous  form  of  Crocco's  Equation  leads  to  an  uncoupled 
expression  for  the  substantial  derivative  of  the  total  enthalpy  in  a  differential  tuiboma- 
chinery  equation  which,  on  integration,  yields  a  novel  form  of  the  'Turbomachinay 
Equation  corrected  for  viscous  non-ideal  flow. 

The  classical  Turbomachinery  Equation  is  an  integrated  expression  arising  from 
Euler's  analysis  of  fluid  torque  and  mechanical  shaft  work  transferred  between  a  fluid 
and  a  rotor  moving  in  that  fluid.  Since  beat  transfer  is  negligible  in  most  turbomachines, 
the  shaft  work  has  been  equated  to  the  specific  total  enthalpy  transfer,  Euler’s  Thrboma- 
chinery  Equation  (s^e  Horlock  [3],  page  77,  Equation  4.3)  as  follows: 


Ah„^A(C/Vu),  0) 

where  Vv  is  the  component  of  the  absolute  velocity,  F,  in  the  direction  U  of  the  rotor 
(or  energy-transferring  device). 

Whether  one  interprets  the  quantity  in  Uk  right  member  of  (1)  as  shaft  woik  or  as  a 
change  in  the  specific  total  enthalpy,  the  right  member  is  a  function  only  of  the  end  points 
of  the  integration.  Therefore,  the  integral  is  an  exact  or  total  integral,  and  it  represents  a 
function  of  state.  However,  in  non-conservative,  dissipating  systems,  attempts  to  express 
the  thermodynamic  process  in  terms  of  the  end  points  alone  leads  to  thermodynamic  in¬ 
consistencies.  An  example  of  such  inconsisteiKies  arises  in  the  consideration  of  a 
propeller  windmilling  on  a  Crictionless  shaft  in  a  moving  fluid.  Since  the  shaft  delivers  no 
work,  a  viscous  correction  is  necessary  whether  the  right  member  of  (1)  represents  work 
or  a  change  in  the  specific  total  enthalpy. 

These  points  are  discussed  in  more  detail  in  terms  of  the  new  differential  equation. 
Also,  an  application  of  the  differential  turbomachinery  equation  is  described  for  a  two-di¬ 
mensional,  ideal,  linear  turbine. 

THE  TRANSFORMA'nON  BETWEEN  ABSOLUTE  AND  MOVING  FRAMES 

In  the  following  discussitm  the  subscripts  v  and  w  represent  the  absolute  and  the 
moving  frame  coordinate  and  vector  values.  (See  Figure  1  and  the  Nomenclature  for  defi¬ 
nitions  of  tpiantities.)  A  frequently  used  relation  connecting  the  absolute  and  moving 
frame  velocities  is 

V^W^U.  (2) 

Following  Spannhake  [4],  spatial  derivatives  in  the  moving  frame  and  time  derivatives 
have  the  following  relationship: 
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Similar  result  are  obtained  for  non-rotating  systems  with  cartesian  coordinates. 


Equations  (3)  and  (4)  define,  in  fact,  the  crypto-steady  criterion  [5],  which,  if  (7  is  ctm- 
stant,  indicate  diat  a  frame  exists  in  which  the  flow  regime  may  be  truly  steady  state. 

Since  the  vector  operator  V  is  independoit  of  time  and  since  vector  operators  are 
independent  of  frame, 

V,»Vw«V.  (5) 

Now  frtnn  Equation  (3)  for  any  static  fruK^tion /, 


where  represents  all  the  position  coordinates  in  the  absolute  frame  and  ^  is  the  mo¬ 
ving-frame  coordinate  in  the  ditecticm  of  motion  of  the  lifting  body  (such  as  the  0 
direction  for  a  rotor).  Equation  (6)  is  an  ^trmsion  of  tihe  usual  crypto-steady  relation. 

lb  appreciate  the  significaiKe  of  the  final  equation  of  this  section,  consido’  an  ob¬ 
server  located  on  a  blade  of  a  rotating  windmill  in  an  infinite  fluid  which  at  infinity 
translates  with  uniform  constant  velocity.  Hie  observer  cannot  detect  any  time-dependent 
changes  in  the  fluid  at  any  given  point  on  die  blade.  However,  on  moving  to  the  absolute 
frame,  the  observer  notes  changes  in  velocity,  pressure  and  temperature  as  each  blade 
passes  (see  Dean  [6]).  Thus,  the  inequality 


indicates  that  in  the  moving  frame  all  partial  time  derivatives  may  vanish  while  the  abso¬ 
lute  partial  time  derivative  may  be  finite. 

Note,  however,  in  contrast  with  the  inequality  (7)  that  the  substantial  derivative  of 
any  static  quantity  at  a  point  is  independent  of  the  coordinate  system  or  frame  and  is 
therefore  an  invariant  addi  respect  to  frame. 
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DERIVATION  OF  THE  DIFFERENTIAL  AND  INTEGRAL  FORMS 


UNCOUPLING  THE  SUBSTANTIAL  TOTAL  ENTHALPY  DERIVATIVE 
IN  THE  MOVING  AND  ABSOLUTE  FRAMES 

Reladonsbips  between  absolute  and  relative  flow  fields  are  given  by  the  Coriolis 
form  of  the  Navier-Stokes  Equations  [3,7-11]. 


—+vvy2-vy  (V  X  r) 

dt 

m^+vwyi-W  y  y  W)+my  W-vuyz 

-Vp  1  „ 

- - —  +  ~V‘3t'  , 

Q  Q 


(8) 


where  Jr'  represents  the  stress  tenscx  excluding  the  pressure  tensor  pd  j  ,  and  U  is  inde¬ 
pendent  of  time.  The  relative  acceleration,  dW/dt ,  is  defined  in  the  relative  frame;  i.e., 
(dW/dt)^. 

Using  the  absolute-fimne  equality  of  (8)  and  the  gradient  form  of  Gibbs’  equation  of 
state,  the  substantial  derivative  of  die  total  enthalpy  is  obtained  in  terms  of  the  partial  de¬ 
rivative  of  the  pressure,  the  substantial  derivative  of  the  entropy  and  the  stress  tensor 

n'  .  SeeWuefal.[7-9]. 


Dho 

Dt 


+  T 


Ds  £ 
Dt  g 


(9) 


At  this  pmnt  non-ideal  analysis  stops  unless  the  partial  derivative  of  die  pressure 
with  respect  to  time  is  resolved  in  terms  of  the  pressure  gradient  [10].  Employing  Equa¬ 
tion  (6)  with  the  static  pressure  as  the  aibitraiy  function  the  result  is 


Dho 

Dt 


U 

Q 


•Vp  +  T 


£i 

Dt 


•(Vw') 


0)  dp 

gde'*' 


r~+-  (v  J?) . 

Dt  Q 


(10) 


Equation  (10)  may  be  simplified  further  using  Wu’s  (reference  [8]  page  91)  defini¬ 
tion  of  the  substantial  derivative  of  the  entropy  in  the  moving  fiame  following  a  particle 
of  fluid  (the  entropy  and  its  substantial  derivative  are  invariant  with  respect  to  frame). 
Thus; 


:  VfF , 
Q 


(11) 


where  q  repments  heat  transfer,  and  the  second  term  is  the  specific  dissipatitm.  Then 
from  (10)  and  (11),  eliminating  die  substantial  derivative  of  the  entropy,  one  obtains  [10] 
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Dho 

Dt 


Q  96  Q  Q 


(12) 


Now  Equation  (8)  may  be  used  to  replace  the  pressure  gradient  in  (12).  From  the  dot 
prodtrct  of  JJ  with  the  moving  frame  equality  in  (8)  the  pressure  gradioit  term  becomes 
(see  Reference  [10] ): 


_  r.  (IT- V)]r+ 2[7. fir  X  IT 

Q  Q  96  9t 

-l7.Vf/V2--l7.(VP). 

Q  (iJ) 

The  gradient  of  the  static  pressure  is  of  course  invariant  in  all  Ihunes,  but  the  mov* 
ing  frame  equatirm  in  (13)  is  convenient  for  later  integration.  Ctnnbining  (12)  with  (13), 
the  substantial  derivative  of  the  total  enthalpy  uncoupled  frmn  the  entropy  and  pressure  is 
obtained. 


-^«[7.i?!+j7.(r.V)r+2l7  firx  r+-V  (^  .«^ +<?  . 

Dt  9t  Q  (14) 

Equation  (14),  the  essential  development  of  this  paper,  is  a  differential  turbomachin 
ery  equation  expressed  in  terms  of  the  moving  frame  velocity  vector  and  q  .  It  is 
universally  applicable  in  any  moving  frame  whether  rotating  or  not. 

Referring  to  Equation  (12),  a  particularly  simple  statement  is  obtained  for  ideal  in- 
viscid  flow: 


Dho/Dt  =  ~U‘  Vp/g  =  - {fodp/dOy/g  (15) 

The  second  equality  in  (15)  is  specific  to  three-dimensional  rotating  machines.  Contrary 
to  Reference  [3],  pp.  7-8,  cmly  the  transverse  jn^essure  gradient  term  survives  in  an  invis- 
cid  system,  indicating  that  neither  axial  nor  radial  pressure  gradients  are  germane  to  the 
calculation  of  specific  total  enthalpy  transfer.  (Mass  flow  rates  are  of  course  a  functiem  of 
axial  or  radial  pressure  gradients.)  the  impulse  stages  of  turbomachines  prove  that  axial 
or  radial  pressure  gradients  need  play  no  role  in  energy  transfer. 

The  ctmclusion  that  there  can  be  no  energy  transfer  without  a  transverse  pressure 
gradient  is  implicit  in  the  logic  of  Euler,  who  insisted  that  only  the  transverse  component 
of  the  momentum  is  effective.  Nevertheless,  the  assumption  that  transverse  pressure  gra¬ 
dients  vanish*is  popular  in  streamline  curvature  methods  of  design  because  it  permits  a 
simplified  two-dimensional  calculation  of  blade  design  [7-9, 11]. 

Equation  (14)  permits  a  revisiting  of  Dean’s  unsteadiness  paradox  [6]  if  one  en¬ 
quires  into  the  possibility  of  rotm  energy  transfer  between  bladeless  discs  or  concentric 


*It  ii  ufued  fiwt  for  an  infinite  number  of  bledea,  dp  becocnec  an  infinitenmal,  and  therefore  dp^  va- 
nidwa  which,  of  coarse,  viotalea  the  rales  of  calcuJnt. 
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cy*indcfs  (t/  =  0) .  (Dean’s  paradox  is  resolved  mosi  simply  by  reference  lo  Equation 
(6)).  A  steady  state  prevails  in  both  the  absolute  and  the  moving  frame,  and  {dWjdt)  va¬ 
nishes.  However,  the  viscous  term  in  (14)  does  not  vanish,  thus  providing  a  mechanism 
for  energy  transfer. 


Dt  Q  Q 


(16) 


Mechanical  energy  transfer  is  accomplished  by  the  second  viscous  work  term  of  the  nght 
monber  in  (16).  It  is  necessary  in  steady-state  systems  that  the  viscous  dissipation  and  the 
heat  loss  balance. 

The  ideal  mechanisms  of  energy  transfo^  contained  within  V  •  Vp/p  and  defmed  by 
(13)  through  (15)  in  their  ideal  limit,  include  a  term  based  upon  the  Coriolis  force  among 
others.  Some  writers  (see  pp.  115-117  of  Reference  [11])  have  suggested  that  the  Coriolis 
term  is  the  necessary  and  sufficient  term  for  energy  transfer.  If  this  were  indeed  true,  then 
axial  machines  and  linear  quasi  two-dimensional  machines  with  vanishing  m,  such  as  the 
wings  of  planes  and  the  sails  of  ships  would  be  useless. 


THE  TOTAL  ROTHALPY  ARGUMENT 

Wu’s  total  relative  rothalpy  relation  [7-9]  (see  also  [12])  with  I  defined  by 

is  given  by; 


Q 


(17) 


Consider  a  system  in  which  the  observer  is  fixed  to  a  blade  of  an  isolated  rotating 
air  screw  or  a  marine  screw  in  an  infinite  uniform  fluid.  No  time-dependence  can  be 
sensed  in  the  moving  frame  and  a  crypto-steady-state  prevails  [6].  The  moving-frame 
partial  derivative  of  the  pressure  with  time  must  vanish.  This  poses  a  great  simplification 
which  justifies  common  practice  [13]  in  the  design  of  marine  screws  where  time  depen¬ 
dence  (and  heat  transfer)  are  generally  ignored. 

Now,  if  viscosity,  heat  transfer  and  the  time-dependent  terms  vanish,  (17)  reduces  to 


WVlmO 

The  well  known  non-trivial  solution  is 


(18) 


V/  -  V(h  +  W^/2  -  U^/2) «  Viko  -  UVa)  -  0  .  (19) 

Equations  (19)  represent  gradioit  forms  of  the  classic  llirbomachinery  Equation  in  the 
moving  and  absolute  frame.  Integratitm  of  the  right  monber  of  (19)  over  a  stream  tube  in 
the  absolute  frame  yields  the  classic  Tlirbomachioery  Equation  (1 ).  Tbe  argument  demtm- 
strates  that  Equation  (1)  is  strictly  true  only  for  isentropic  and  steady-state  flow  in  a 
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stream  tube.  However,  in  systems  with  dissipation,  it  is  necessary  to  account  for  non-ideal 
effects  in  (1). 


MOVING  FRAME  TIME  DEPENDENCE 

Consido'  in  Equation  (14)  a  small  sinusoidal  moving-frame  oscillation  imposed  on 
the  relative  velocity  W  leading  to  a  new  velocity,  W  .  The  ratio  of  the  magnitude 
of  time-^pendent  W'  to  time-independent  W  will  be  arbitrarily  fixed  at  ^  0.05  .  The 
imposed  frequency  of  the  time-dependent  IT*  is  at  least  an  order  of  magnitude  greater 
than  the  blade  velocity. 

'Wldi  Reynolds  averaging,  the  partial  derivative  of  W  with  respect  to  time  will  van¬ 
ish.  However,  the  Reynolds  average  over  non-linear  terms  will  lead  to  non-vanishing 
Reynolds  stress  terms  which  will  not  be  addressed  here.  As  a  result  of  the  Reynolds  time¬ 
averaging  process.  Equation  (14)  may  be  written 

-l7-((ir- V)  W\-¥2V-Sf  X  (R^)  +  (l/p  V-Jt'  irj  +  l  q  )  ,  (20) 

where  the  braces  represent  the  Reynolds  time  average  and  W  and  tt*  include  the  time- 
dependent  W  .  This  notatirai  will  not  be  used  again. 

In  the  Re}molds  time-averaging  process,  the  acceleration  term  drops  out  as  ex¬ 
pected.  Elimination  of  the  acceleration  term  simplifies  the  integration  of  (14).  As 
mentioned,  heat  transfer  is  generally  not  significant  in  most  rotors  [11,13],  and  in  partic¬ 
ular,  those  associated  with  pumps  and  marine  screws,  and  it  will  be  ignored  here. 

Now,  note  that  arbitrary  samples  of  turbulent  flow  over  short  (relative  to  the  blade) 
periods  may  be  written  as  a  Fourier  summation  of  regular  sinusoidal  oscillations.  For  the 
limited  picture  developed  here.  Equation  (20)  is  again  applicable  according  to  the  above 
argument 


INTEGRATION  OF  THE  TOTAL  ENTHALPY  RATE 

A  proper  test  of  (14)  would  be  whether,  on  integration  over  the  rotor  blade-to-blade 
flow,  it  would  predict  a  total  enthalpy  transfer  compatible  with  that  of  the  classic  Thrbo- 
machinery  Equation  (1).  Therefore,  the  integration  of  (14)  for  crypto-steady  conditions 
will  be  p^ormed  as  a  test  in  the  three-dimensional  domain.  Then  an  application  of  the 
new  equation  will  be  developed  in  a  two-dimensional  linear  turbine. 

In  the  integration  process  it  will  be  assumed  that  the  flow  may  be  divided  into 
streams  which  pass  between  a  given  pair  of  blades.  In  the  rotating  frame  the  streamtube 
walls  are  fixed  steady-state  walls  associated  with  a  steady-state  mass  flow  rate  m  which 
may  consist  of  radial  and  axial  mass  flow  components. 

DERIVAnON  OF  THE  INTEGRAL  FORM  FROM  THE  DIFFERENTIAL  FORM 

The  substantial  derivative  of  the  total  enthalpy  represents  power  transfer  at  any 
point  in  a  flow  field  between  the  fluid  and  a  body  moving  at  constant  speed  in  that  fluid. 
To  obtain  the  total  enthalpy  transfer  rate  over  the  entire  field,  it  is  necessary  to  integrate 
(14)  thus: 
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III'? 


Q  [Z7-(^  V)^+2l7-fir  X  W]  rdrdBdz 


+  ^  j  jv  W)  rdrdBdz  , 

where  the  tecoai  term  of  the  right  member  is  the  viscous  term.  The  first  term  of  the  right 
member  of  (21)  is  hie  tangential  component  of  the  convecthm  term  (^Hained  on  dot  multi¬ 
plication  widi 

qU-(W''V)W 


\  dr  r 


9Wg  ^  W,Ws\ 
dO  *dz  r  / 


(22) 


The  first  and  fourth  terms  of  (22)  will  be  combined  in  an  integral  identified  by  /|,4  thus: 

arWfl 


liA 


mWr  9rV 


■rdrdBdz 


(23) 


where  we  have  used  the  mean  value  theorem  to  take  W9  (r)  outside  hie  double  integral. 
Thus: 


^e(r) 


W#  dBdz  . 


(24) 


The  factor  in  parenhieses  in  the  right  member  of  (23)  is  the  radial  mass  flow  at  any  point 
r 


mf(r) 


QW^d9dz~mfAr) 


where  Mr)  is  the  radial  firachon  of  the  mass  flow  rate  m.  Then 


(23) 
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(26) 


Now,  io  the  steady-state  the  mean  value  of  /<r)  ,fr  ,  is  a  constant  given  by 

^ejch 

fr  l  MU^0)diuWa)/MU^a)  ,  <27) 

and, 

/M-i<A<W-iM(tWs)r,-(CAV,)^l  (28) 

where  ihr  is  the  mean  value  of  .  Using  die  same  arguments  as  used  for  the  radial 
mass  flow,  the  integral  of  the  third  term  of  the  right  member  of  (22)  may  be  written  to 

show  the  axial  mass  flow  rate  m,(2)  exjdkitly.  Widi  the  velocities,  W, ,  Ws ,  and  0  av¬ 
eraged  over  r  and  $ ,  one  obtains  the  steady-sude  mean  value,mx  ,  and 


h  -  .  (29) 

The  second  term  of  the  right  member  of  (22)  provides  an  integral,  /2 ,  which  con¬ 
tains  the  tangential  kinetic  energy. 

'‘‘I 

Since  the  tangential  velocities  at  the  blade  walls  are  the  blade  velocity,  the  integral  va¬ 
nishes. 

Now  identifying  the  second  term  of  the  first  integral  in  the  ri^t  member  of  (21  as 
Is,  we  may  write  using  die  above  argummits 

/s-  j  I J 2ewhW^drd9dz^mA(l/\-((/^)r,]  .  OD 

Summing  die  components  of  integratkm  (28)  throu^  (31) 
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(32) 


+  m,A,(^«)  +  (V/5)  , 

where  Ar  sod  A|  refnesent  die  change  along  r  and  z  respectively,  and  (VIS)  is  the  inte¬ 
grated  viscous  tenn. 

Now  adding  which  is  zero  to  (32), 

I  j  |p^dr-inA[(t^+^^a)£^] 

+ mA[(i^ + +  (V/5)  -  .  (33) 

In  (33)  the  terms  W0  are  averaged  over  6  and  z  in  the  first  term  and  over  r  and  in  the 
second  term.  Finally, 

Aho  ~fAA09)  +/A(t/^a)  +  (V/S)/m  -  A(t/^s)  +  (V/S)/m  .  (34) 

In  the  right  member  of  (34)  die  term  (V7S)  contains  die  viscous  work  as  well  as  the 
dissipation,  which  may  provide  a  net  outflow  of  eiMrgy  in  viscously-coupled  biadeless 
devices.  The  delta  term  which  approximates  the  right  member  of  Euler's  Equation  (1 ) 
may  be  ctmsidered  a  pseudo-ideal  term. 

Equation  (34)  exhibits  srane  similarity  with  the  ideal  Equatkm  (1).  The  fact  that 
mathematical  ^tial  averaging  processes  have  been  enqiloyed  suggests  that  the  pseudo- 
ideal  term  cannot  be  a  function  of  the  end  poims  alone.  Therefore,  the  pseudo-ideal  term 
cannot  be  ideal  unless  the  system  is  ideal.  Similar  conclusions  apply  to  die  total  relative 
rolbalpy  [12] ,  which  depends  upon  the  pseudo-ideal  term. 

It  is  easy  to  combine  the  axial  and  radial  flow  terms  in  (34)  by  letting  the  operator 
A  vary  in  both  r  and  z .  Hie  derivation  siqtpoits  the  conchision  that  Equation  (14)  is 
indeed  a  differential  turbomachmery  equation. 

A  two-dimensional  application  and  test  of  die  differential  form  (14)  on  an  ideal  lin¬ 
ear  device  where  the  solution  is  known  preci^y  will  now  be  examined. 

THE  SUBSTANTIAL  TOTAL  ENTHALPY  RATE  IN  A 
TWO-DIMENSIONAL  DEVICE 

An  infinite  circular  cylinder  widi  bound  circulation,  as  shown  in  nguie  2,  is  an  ele¬ 
mental  linear  turbine.  It  may  be  considered  as  an  infinite  sail  on  a  sailboat  or  an  infinite 
wing  on  a  sailplane.  The  device  extracts  energy  from  the  ideal  inviscid  working  fluid. 
Work  is  performed  on  the  sailplane  (fibred  to  a  vertical  rail)  by  rr  Ising  its  height  at  uni¬ 
form  sp(^  V  agaiiat  gravity.  Work  on  the  sailboat  is  performed  by  moving  die  boat  at 
uniform  speed  U  adiich  elevates  a  weight  attached  at  minus  infinity  by  an  infinite  tether. 
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In  the  moving  frame  the  apparent  velocity  of  the  ideal  working  fluid  at  infinite  distance  is 
Wo .  The  relationship  between  the  absolute  and  moving  coordinate  system  and  the  velo¬ 
cities  is  given  by  the  transformation  of  Figure  1. 

Preston  [14]  computed  energy  transfer  in  two-dimensional  arrays  of  ideal  vortex 
points  by  transforming  their  potential  functions  to  the  absolute  frame  where  they  exhibit 
time  dependence.  The  study  involves  circular  cylinders  of  finite  dimensions  with  bound 
vorticity  which  can  be  transfonned  into  two-dimensional  airfoils.  Energy  transfer  is  de¬ 
veloped  in  steady-state  moving  frames  by  integrating  Equatitm  (14)  over  all  space. 

THE  TOTAL  ENTHALPY  TRANSFER  RATE  BASED  ON  AERODYNAMICS 

Since  the  flow  field  is  ideal,  the  flow  domain  may  be  described  by  a  potential  func¬ 
tion  m  its  conjugate  stream  function.  The  lift  is  therefore  the  ideal  lifting  force.  L,  of  the 
Kutta-Joukowsld  Equation  given  by 


L^qWJ  ,  (35) 

where  F  is  the  scalar  circulation.  The  units  are  force  per  unit  length  of  cylinder.  In  the 
absolute  and  moving  frame  the  lift  component  Ly  directed  parallel  to  the  y  axis  of  Figure 
2  is  given  by 


Ly  ■  QWo,^  «  ,  (36) 

where  the  subscript  x  represents  the  x  component.  Recalling  that  is  the  velocity  of  mo¬ 
tion  of  the  device  (sail  or  wing  or  rotating  cylinder)  as  perceived  in  the  absolute  frame, 
the  power  is  the  product  of  U  and  Ly. 


Power  I  unit  length  ■  pl/lPoxF  .  (37) 

Since  we  assume  that  there  is  no  beat  rate. 


DHq’ 

Dt 


-QUWoj:  , 


(38) 


where  Ho*  is  the  total  mitfaa^iy  of  the  system  per  unit  length  d  lifting  surface.  Equation 
(38)  is  the  anticipated  relationship  which  should  ultimately  be  developed  from  the  differ¬ 
ential  form  (14). 


THE  STREAM  FUNCTION,  VELOCTTY  AND  RELATIVE  ENTHALPY 
IN  THE  FRAME  OF  THE  BLADE 

Since  ideal  flow  has  bemi  uisumed  in  the  moving  frame  of  the  blade,  the  stream 
function,  V' ,  is  the  usual  function  modified  for  motion  along  the  y  axis. 
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v>  *  ~  Woyil  -  a^/r^)r  cos  6  +  WpA}  -  a^fr^y  sin  6 


+  (r/2Jt)ln(r/a)  .  (39) 

The  constant  a  is  the  radius  of  the  cylinder.  The  cartesian  velocity  components  are 
obtained  by  the  usual  transfonnaticm  [10]  as  follows: 

r  y 

2jr  (?+7)  ’  (40) 

.  <T2^^(x2-y2) 

__L 

2w  (?T?)  *  (41) 

Now  the  relative  vorticity  must  vanish  because  potential  flow  cannot  have  vorticity. 
A  check  of  the  vorticity  in  the  relative  frame  shows  tltot  indeed  it  vanishes.  Also,  the  ti* 
me<dependent  term  vanishes. 

THE  SUBSTANTIAL  TOTAL  ENTHALPY  DERIVAnVE 

In  the  linear  two<diinensioiud  system,  tte  differential  form  of  the  turbomachinery 
equation  (14)  is  simplified  because  the  rotation  vanishes. 


Since  the  vorticity  vanishes 


w-^^vwyi , 


Dt 


UdW^ 

2  dy*  ■ 


(44) 


The  integrated  substantial  mtal  enttia^  rate  per  unit  length  [where  die  subscript  cm  y  in 
(44)has  beendrt^iped]  is 


Integration  of  (45)  wiU  be  performed  over  an  space  per  unit  length  z  of  dm  blade.  The 
choice  of  time  is  immaterial  since  the  fluid  dynunics  are  steady  state  in  die  moving  frame 
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and  the  thermodynamic  rates  over  all  space  are  invariant  with  time.  It  is  understood  that 
the  integration  applies  only  to  the  fluid  domain  and  that  boundaries  at  solid  walls  are  ob¬ 
served. 

Now  without  going  into  details  [10]  the  integral  of  the  total  derivative  in  (45)  is 


--i:( 


AUW„  x{a^  - x2)>/2(2x2  _  ^2) 


+ - ^ — 2 - + - r - 

jta^  a* 


y 


Note  that  only  odd  terms  in  y  make  any  contribution  to  (46).  Since  the  first  and  fifth 
terms  cancel,  only  four  terms  remain.  The  integration  with  respect  to  x  is  performed 
through  a  transformation  employing 


jr  =  flcos0  , 


with  integration  limits  given  by 


9  •‘it  when  ar*-a  , 
^  *  0  when  x  «=  a  . 


Making  the  substitutions 


co^Osin^SdG 


-%WaxUa  I  cosdsin^^dO 


f”.i 

»  J. 


sin^ddi9 


•5^1  8in2»(l-2cos^«)<»  j  . 
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In  (50)  the  second  integral  makes  no  contribution  because  it  is  antisymmetric.  Tbe  first 
integral  cancels  the  s^ond  term  in  the  last  integral  to  yield  from  the  surviving  terms 


DH  * 

Equation  (SO)  is  identical  with  (38)  and  this  result  illustrates  a  useful  ai^lication  of 
the  diffemitial  form  and  constitutes  confirmatitHi  of  tbe  validity  of  the  differential  turbo- 
machinery  Equation  (14).  Pdr  the  linear  case,  the  energy  transfer  rate  of  tbe  rotor  is 
proportional  to  the  component  of  the  kinetic  energy  gradientparallel  to  the  moving  rotor. 
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